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ABSTRACT 



u: 

^ . We set out to determine the ratio, giR, of rest-frame 8-1000-/im flux. Sir, to monochromatic radio flux, 5i.4ghz, for galaxies 
5— ( ' selected at far-infrared (-IR) and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust 
^ temperature, T^, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation 
, (FIRRC) is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-/im and 
' ' . 1.4-GHz samples, obtained using Herschel and the Very Large Array (VLA) in GOODS-North (-N). We determine bolometric IR 
' output using ten bands spanning Aobs = 24 — 1250 /^m, exploiting data from PACS and SPIRE (PEP; HerMES), as well as Spitzer, 
SCUBA, AzTEC and MAMBO. We also explore the properties of an LiR-matched sample, designed to reveal evolution of giR with 
redshift, spanning log Lir = 11-12 L© and z; = — 2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies 
in GOODS-N, we see tentative evidence of a break in the flux ratio, giR, at 11.4 ~ 10^^'^ WHz~^, where active galactic nuclei 
(AGN) are starting to dominate the radio power density, and of weaker correlations with redshift and Td. From our 250-/im-selected 
sample we identify a small number of far-IR-bright outliers, and see trends of giR with 11.4, Lir, Td and redshift, noting that some 
of these are inter-related. For our LiR-matched sample, there is no evidence that qir changes significantly as we move back into 
the epoch of galaxy formation: we find giR oc (1 -I- z)'' , where 7 = —0.04 ± 0.03 at 2; = — 2; however, discounting the least reliable 
data at 2 < 0.5 we find 7 = —0.26 ± 0.07, modest evolution which may be related to the radio background seen by ARCADE 2, 
' perhaps driven by <10-pJy radio activity amongst ordinary star-forming galaxies at 2; > 1. 

Key words. gala:xies: evolution - galaxies: starburst - infrared: gala:xies - submillimeter: galaxies - radio continuum: galaxies 
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1. Introduction 119891: iLisenfeld et al.lll996l ). dust absorbs all of the ultra- 

■ violet radiation from massive stars, re-radiating this en- 

■ - • For samples of local galaxies - on galactic and -100- ergy in the far-IR, and when those massive stars explode 

pc scales - there is a good correlation between far-IR as SNe they generate cosmic-ray electrons which lose all 

and radio emission (dc Jong ct al .' [19851; iHelou et al.||l985|; then energy in the radio regime, mainly via synchrotron 

[Condon ct al. 1991; Yun ct al. 200lJ). The correlation spans emission. A balance is thereby achieved between far-IR and 

many orders of magnitude in luminosity, gas surface den- ^adio emission, assuming that the starburst timescale is suf- 

sity and photon, cosmic-ray and magnetic energy density, flciently long (>10^yr) 

and arises because the far-IR and radio wavelength regimes Traditionally, and L^adio are both employed to de- 
share a common link with luminous, massive stars and then termine star-formation rates, and the far-IR/radio flux den- 
end products - dust, supernovae (SNe) and cosmic rays gj^y j-atio has been useful when estimating the redshift or 
In the simplest models (dubbed 'calorimetry' - e.g. |VoelU Th of a distant starburst, or when defining samples of AGN 
(Condon 1992: Carilh k Yun '1999; Ivison et al. '200^ 111 



*Herschel is an ESA space observatory with science instru- 2003; Chapman ct al. 2005; Donley ct al. 2005) or probing 
ments provided by European- led Principal Investigator consor- magnetic field strength (jThompson et all 120061). For these 
tia and with important participation from NASA. reasons, and because of recent observational advances at 
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both far-IR and radio wavelengths, there has been a deluge 
of FIRRC-related work recently, exploring why the correla- 
tion exists and whether it contin ues to hold at progressively 



larger look-back t imes ([Garret t 2002; Appl eton et alJ 



Ibar et al. 


20081: Sevmour et al.. 


Sargent et al. 2010a). Prevailing 


20101) sugj 


jests that variations in 1 



2004 : 



201C; 



should be virtually unavoidable and that the FIRRC thus 
arises due to a mysterious combination of effects involving 
bremsstrahlung, inverse Compton cooling, ionisation and 
the relative fractions of primary/secondary cosmic-ray elec- 
trons/protons, as well as the critical sy nchrotron fre quency. 

Aside from the modelling work of iLacki et al.l . recent 
advances in this field have included the use of luminosity- 
matched samples (between high and l ow redshift) to better 
probe evolution with look-back time (jSargent et al.l[2010b[ ) 
and the use of measurements spanning the far-IR and radio 
wavebands to avo id assumpt ions relating t o k co rrections 
(jlvison et al.ll2010[ ). although [Calzetti et al.l (|2010[ ) have ar- 
gued that bands beyond 24 fim contain a contribution from 
dust heated by stars from previous episodes of star forma- 
tion and so we might not necessarily expect the correlation 
to improve. In this paper we introduce flux-limited 250-/im- 
and 1.4-GHz-selected samples of galaxies from Herschel and 
the VLA, as well as a luminosity-matched sample selected 
at 24 /im, spanning z = — 2, and determine their spec- 
tral energy distributions (SEDs) spanning the entire far-IR 
spectral region. We then investigate the FIRRC from the 
perspectives of the 24-, 250- //m- and radio-selected samples. 

2. Sample selection and data analysis 

In this p aper we present resu lts from observations with 
Herschel (jPilbrattl let al. 2010D . The SPIRE instrument, 
its in-orbit p erformance, and its scientific capabilities are 
described by iGrifhn et all (|2010D . and the SPIRE astro- 
nomical calibra t ion m ethods and accuracy ar e outlined in 
ISwinvard et"all (|2010D . PACS is described bv iPoglitschI et 
al. (2010). 

Our datasets are drawn from the common area observed 
by PACS and SPIRE at 100, 160, 250, 350 and 500 /xm as 
part of HerME^ (Oliver et al., in preparation) and PEP 
(Lutz et al., in preparation) in the GOODS-N field, prior 
to acquisition of data for GOODS-Herschel. GOODS-N has 
a lso been observed with the VLA at 1.4 GH z (1.7" fwhm 
- iBiggs fc Ivisonll2006l: [Morrison et al.ll2010D and Spitzer at 
24, 70 and 160 /xm; we make use of th ese data, as well a s 
the 850-, 1100- and 1250-jUm images of 'Bor vs etall (|2003D . 
iPerera et all (|2008D and[Grcvc et al. (2008). 

We employ three GOODS-N galaxy samples, aU selected 
above a signal-to-noise threshold of 5 a: 

1. 128 galaxies selected at 250 /im, without priors, with 
<5'250/im ^ 20mJy (Fig.[TJ Smith et al., in preparation); 

2. 247 galaxies selected at 1.4 GHz (Fig. [T|) with a 
'5'i.4GH z limit of ~20/xJy, 137 with spectroscopic red- 
shifts ( Barger et al.l [20081) . the remainder with photo- 
metric redshifts {{z) = 0.94; interquartile z, 0.56-1.76); 

3. a LiR-matched sample of 652 so urces spanning z = 
- 2, sel ected initially at 24 ^m ([MagneUi et al.ll2009t 
iBerta et'al . 2010) then filtered to cover only the decade 
of LiR between 10^^ lO^^LfT^ (LIRGs), where Lir i s 
determined using the models of lCharv &: Elbad (|200lD . 

^ hermes.sussex.ac.uk 
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Fig. 1. 250-Aim image of GOODS-N, with >5-a 250-fim [ra- 
dio] sources marked by circles [crosses]. The PACS region is 
also indicated. Here, a includes instrumental and confusion 
noise combined in quadrature, where CTconf /cinstr 5 (cf. 
~2 for BLAST). 



Far-IR and submm flux densities for the three samples 
are determined using images convolved with appropriate 
point spread functions. Sm is calculated by integrating un- 
der the well-sampled SEDs. Monte-Carlo simulations are 
used to assess the uncertainty in S'ir. The formal error 
on S'24pm was boosted by 3x to account for the uncertain 
shape of the SED between rest-frame 8-70 /im. A modified 
blackbody fit to the measurements beyond 24 /im (with the 
emissivity index, /3 — 1.5) was used to determine T^. 

For sample (1), additional procedures are implemented 
to define a clean sample, free fr om blends: following the pro- 
cedure of iDownes etaLl (|1986[) . 107/128 sources are found 
to have secure (P < 0.05) radio identifications (ids) within 
a search radius, r = 10"; we discard the remainder. To 
avoid using those sources most severely affected by blend- 
ing, we further discard those with more than one radio emit- 
ter within r, leaving 65 sources. Of the galaxies without a 
secure radio id, three have no plausible radio ids within 
r: a potentially interesting sub-sample. Measurements are 
made at the radio positions for the 65 sources with secure, 
unambiguous ids, and at the 250- /im positions for the three 
sources without radio emission. 

For sample (2), far-IR and submm measurements are 
made at the radio positions. 

For the LiR-matched galaxies (sample 3), median stack- 
ing is used to me asure Sm. and 5j^4GHz : we follow the proce- 
dure outlined by llvison et al.l (|2010t) . Fluxes are calculated 
from 31^-pixel^ stacked images in the ten available filters 
and S'ir is determined as before. 



3. Results and conclusions 

qiB. as utilised here is the logarithmic ratio of the 

rest-frame 8-1000-/i,m fiux, 5ir, and the 1.4-GHz flux 

density, 5'i.4ghz, such that qm = logio [(S'ir/3.75 x 
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Fig. 2. giR versus Sin for those 250-/xm-selected galaxies 
(sample 1) with secure, unambiguous radio ids. Those with- 
out plausible radio ids are plotted as stars. The dashed line 
is the median, giR = 2.40; the shaded region represents 
±2ag (o-g = 0.24). 

Table 1. Trends. 

qm trend Spear- Signifi- 

man p -cance 
Sample 1 (250-/im-selected galaxies with redshifts): 
(5.02±0.18) - (0.105±0.008)logll.4 -0.48 99.8% 
(6.09±0.33) - (0.092±0.008)logLiR -0.32 95.6% 
(2.61±0.02) ~ (0.081±0.007) (l-f^) -0.33 96.0% 

(2.76±0.03) - (0.008±0.001) T,i -0.33 96.1% 

Sample 2 (radio-selected galaxies with redshifts): 
(4.92±0.21) - (0.101±0.009)logll.4 -0.27 99.9% 
(2.74±0.35) - (0.007±0.009)logLiR +0.07 69.1% 
(2.55±0.02) - (0.047±0.010) (l+z) -0.15 96.6% 
(2.60±0.02) ~ (0.002±0.001) Td -0.16 89.2% 



10i2Wm-2)/(5i.4GHz/Wm-2Hz-i)], where 5i.4GHz is k- 
corrected assuming Sv ocu", with a = —0.8. 

We begin with sample (1), those selected at 250 //m: giR 
is not a strong function of S'ir (Fig. [5]), nor of S'i.4ghz- We 
see no evidence of contami nation by r a dio-lo ud AGN, con- 
sistent with the findings of lYun et al.l (|200lh . Some galax- 
ies stand out as potentially far-IR-bright: these include the 
three galaxies without plausible radio ids, two of which are 
detected at 70 and/or 160 /im, so are likely at low redshift 
with their radio emission resolved away. 

Only 39/65 sources with unambiguous radio ids 
have redshifts (20 photometric, 19 spectroscopic; {z) = 
0.98; interquartile z = 0.46-1.52, similar to sample 2). 
Nevertheless, this sub-sample allows us to explore correla- 
tions between and luminosity, redshift and Td- We find 
significant (>95% confidence ~ Table [T]) trends for lower 
giR, amongst the most radio- and far-IR-luminous galaxies, 
and the warmest and most distant, though these parame- 
ters are likely inter- related. The dependence of giR on 11.4 is 
the strongest and likely reflects the influence of low-radio- 
power AGN, of which more later; that of giR on Ljr is 
more puzzling, perhaps reflectin g the dependen ce of Ljr 
on redshift and/or fe.g. , Chapman et al.ll2005h . or selec- 
tion effects (since this trend is not seen for sample 2 - see 
Table [I]). 
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Fig. 3. giR versus redshift for our radio-selected galaxies 
(sample 2), in five bins of X-corrected 11.4, plus the full 
sample. Values of giR for sample 3 are shown (as circles) 
for comparison. 



Fig. [3] shows giR versus redshift for our radio-selected 
galaxies (sample 2), split into five log-spaced bins of 11.4. 
Does giR evolve with redshift? One might conclude that 
it does, based on the bottom panel of Fig. |3l where giR 
oc (1 + z)T, with 7 = -0.05 ± 0.01 (Table [1]). However, 
we must be aware of some strong selection effects which 
make this evidence unreliable: radio emission can be due to 
an AGN and several radio-loud objects with low values of 
giR are obvious in Fig. |31 Such AGN are more common at 
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Fig. 4. Median giR versus 11.4. The local luminosity func- 
tions of starbursts and AGN are shown (jMauch fc Sadlen 
[2001 . 




formation at z ~ 2, with 7 ~ —0.04 ± 0. 03 where qm 
( X (1 + zy, consistent with the findings of I Sargent et al.l 
()2010bD . If we discount the z < 0.5 data, which comprise 
only 16 galaxies which are not well matched in Lm to the 
higher redshift bins, we find 7 = — 0.26±0.07. This is si milar 
to the 7 = -0.15 ± 0.03 found bv llvison et al] (|2Q10D who 
noted reports that evolution in giR could b e related to the 
radio background s een by ARCADE 2 (Fi xsen et al.|[2010l: 
ISeiflFert et al.ll2010l ). Our sample, with (S'i.4ghz) ^ 10 /iJy 
at z > 1, is consistent with the idea that evolution of the 
FIRRC might be driven by <10-^Jy ra dio activity amongs t 
ordinary star- forming galaxies at z > 1 (ISingal et al.ll2010l ) . 
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Fig. 5. Liji (dots; left axis) and giR, (red ci rcles; right axis) 
- the former determined via the models of ICharv fc Elbaa 
(|200l|) - versus redshift for our LiR-matched sample. 
The luminosity bounds and redshift bins (dashed lines), 
the number of galaxies in each bin and their predicted 
(measured) (logLjR) and measured (^ir) are all shown. 
The shaded are a rep resents a ±la p rediction for giR 
(jSwinbank et al.ll20od: llvison et al.ll2010[ ). 

z ~ 2 than today (e.g. IWall et aLll2005[ ): moreover, radio 
emi ssion from faint starbursts (with a = —0.8, although 
see llbar et aTl |2010[ ) becomes more difficult to detect at 
higher redshifts, such that the fraction of radio-loud AGN 
in a fiux-limited sample will rise, driving down giR. Indeed, 
Fig. U shows tentative evidence of a break in {qm) at 11.4 
10^^-^ WHz^-'^. One might also expect radio-loud objects 
(those with low (/ir) to contain warmer, AGN-heated dust, 
giving rise to the weak trend (89.2% confidence - Table [1]) 
of decreasing giR with increasing Tj. 

Finally, we turn to our LiR-matched galaxies (sample 
3), illustrated in Fig. [S] The Sz = 0.5 bins provide sig- 
nificant numbers of objects at near-constant Lir spanning 
z = — 2. As well as being matched in Lir, there is another 
key difference betw een our new sample and that used by 
llvison et al.l (|2010[ ): although the new sample is based ini- 
tially on a flux-limited 24-/im catalogue, the final selection 
is based on Lir, with model-dependent extrapolations from 
the mid-IR (accurate to < 2x across all bins - Fig. [5]). This 
should lead to less contamination by AGN at the blue end 
of the rest-frame 8-1000-/im band, where the relative con- 
tribu tion to SiB, can be substantial (Figure 11 - llvison et al.l 
[20T0I) . Using our new sample, there is no strong evidence 
that giR changes as we move back into the epoch of galaxy 
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